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Surface plasmon polaritons (SPPs) are coherent electromagnetic surface waves trapped on an 
insulator-conductor interface. The SPPs decay exponentially along the propagation due to con¬ 
ductor losses, restricting the SPPs propagation length to few microns. Gain materials can be used 
to counterbalance the aforementioned losses. We provide an exact expression for the gain, in terms 
of the optical properties of the interface, for which the losses are eliminated. In addition, we show 
that systems characterized by lossless SPP propagation are related to PF symmetric systems. Fur¬ 
thermore, we derive an analytical critical value of the gain describing a phase transition between 
lossless and prohibited SPPs propagation. The regime of the aforementioned propagation can be 
directed by the optical properties of the system under scrutiny. Finally, we perform COMSOL 
simulations verifying the theoretical findings. 
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I. INTRODUCTION 

A light-matter interaction called surface plasmon po¬ 
laritons (SPPs) has gained the scientists’ interest due 
to its unique properties, such as control of electromag¬ 
netic energy in subwavelength scales high sensi¬ 
tivity in dielectric properties , negative refraction 

and hyperbolic wave front 0,0]. SPPs have been applied 
in nanophotonics, imaging, optical holography, nano an- 
tennas, biosensing, integrated circuits and metamaterials 

IMt 13 ■ Important progress has been made in plasmon- 
ics with two-dimensional materials, such as graphene and 
black phosphorus, where the plasmonic properties can be 
tuned by u sing chemical doping or applying external gate 
voltage [3, [l3l - fl8l| . Moreover, plasmonic lenses, waveg¬ 
uides and meta-materials based on graphene have already 
been applied [ 1 , [l^[2^. Last but not least, multi-layers 
structures have been created by stacking two-dimensional 
crystals one on top of another providing surprising elec¬ 
tronic and optical features [13, [13, [HI . 

Near plasma frequency ujp, the electrons on the sur¬ 
face of metals or semiconductors are free to move sus¬ 
taining collective oscillations [i, [lii^. The coupling 
between light and electron oscillations allows the creation 
of Transverse Magnetic (TM) Electromagnetic (EM) sur¬ 
face waves, namely surface plasmon polaritons (SSPs). 
From the mathematical point of view, SPPs are surface 
waves bounded along the interface between two mate¬ 
rials with sign reversed dielectric permittivities, i.e. a 
dielectric-conductor interface, and their EM field decays 
exponentially away from the interface (evanescent waves) 

[ 3 , Mill- 

In this work, we focus on plasmonic waveguides formed 
by a planar interface which consists of two semi-infinite 
layers with reversed sign permittivities, namely a dielec¬ 
tric and a metal. The dispersion relation which char¬ 


acterizes the SPPs propagation can be determined by 
Maxwell Equations (ME) as [ 3 , 0 ) 113113 | 


P = ko n = ko 






( 1 ) 


where fco = w/c is the free space wave number of the 
incident excitation light of angular frequency w, n is the 
plasmon effective refractive index, Sd and Sm the permit¬ 
tivity of dielectric and metal, respectively, and c is the 
speed of light in vacuum. 

SPPs decay exponentially along the interface as well, 
due to the metal losses. In mathematical language, the 
metal losses are described by a negative imaginary part 
in the permittivity function of the metal, i.e. Em = 
—Sim — i£ 2 m, where £imF 2 m > 0. Consequently, the 
SPPs wave number /3 becomes complex, viz. /3 = /3'-|~i/3", 
where the imaginary part accounts for losses of SPPs en¬ 
ergy. The imaginary part 5[/3] of Eq. ([T]) determines 
the characteristic propagation length L, which shows the 
rate of change of the e nergy a ttenuation of SPPs along 
the propagation axis 0 , l23l 1 ^ , that is 


2^PY 


( 2 ) 


Gain materials, rather than passive dielectrics, have 
been used to reduce the losses in SPPs propagation. 
These active materials are characterized by a complex 
permittivity function, i.e. Ed = Eid + iE 2 d with Eid, E 2 d > 
0, where the imaginary part acconnts for gain, that is, 
the dielectrics give ener gy t o the system counterbalanc¬ 
ing the metal losses [1 in, [1313. In addition, active 
dielectrics have been used for exploring VF symmetry 
in optical systems [33 - 1^ characterized by the condition 
that n{—x) = n*{x), where n and n* the refractive in¬ 
dex and its complex conjugate, respectively; x denotes 
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the spatial coordinate along the interface. Metamateri¬ 
als with V'T symmetric effective refractive index can be 
constructed the combination of gain dielectrics and 
loss metals [3ll - [^ . What makes VT symmetric media 
interesting is that they allow control over EM field by 
timing the gain and loss of the materials. 

It has been already demostrated in [3, [H, [H, 113 that 
for a certain value of gain, the losses in SPPs propagation 
may vanish. Consequently, the SPPs propagation con¬ 
stant (3 as well as the effective refractive index n become 
real and therefore the VT symmetry is satisfied, since n 
does not exhibit any spatial dependence along the inter¬ 
face. Furthermore, Eq. o states that a VT symmetric 
n leads to infinite propagation length, viz. lossless SPPs 
propagation. 

In the present work, we investigate theoretically and 
numerically the VT symmetry in active plasmonic sys¬ 
tems. In Section m we provide an explicit expression 
of the gain, namely epT, for which the losses in SPPs 
propagation have been eliminated. In addition, we find 
a critical value Sc of e-p-p, where SPPs wave number P 
and the effective SPPs refractive index n shift from real 
to imaginary regime, subsequently the Sc is a VT sym¬ 
metry breaking point. It is remarkable that it is a steep 
phase transition from lossless to prohibited SPP propa¬ 
gation, which offers the opportunity to control whether 
SPPs propagate or not by tuning the optical properties 
of the interface. In Section imi we apply the theoreti¬ 
cal results derived in the previous Section on interfaces 
comprised of active dielectrics and Drude metals. In Sec¬ 
tion IlYl we proceed with numerical simulations by solv¬ 
ing the full system of ME in the frequency domain by 
using the commercial multiphysics software COMSOL, 
and we show that lossless SPPs propagation correspond¬ 
ing to VT symmetry can be achieved in the presence of 
gain dielectrics. Finally, concluding remarks are offered 
in Section |3 


II. VT AND CRITICAL GAIN 


In this Section, we calculate the exact expression of 
the dielectric permittivity gain counterpart e2d, for which 
the SPPs propagate without losses in the dielectric-metal 
interface. 

Plugging the complex structure of the dielectric and 
metal permittivity into Eq. o, function n can be written 
in the ordinary complex form as [33| 

+ i sgn(y) ^ (3) 


n = 


y/x'^ +y'^ + X 


where sgn(?/) is the discontinuous signum function and 


^ld\\£r, 


^Im ll^dll 


y ■■= 


ll^d + £m IP 

-£2m||ed||' 


£2d||£r; 


IPd + £m IP 


(4a) 

(4b) 


with ||e,|| denoting the norm of the complex number e*. 

Considering the plasmon effective index n in Eq. © 
in the (x, j/)-plane, we observe that a lossless SPP prop¬ 
agation, i.e, 3[/3] = 0, is warranted when the conditions 
2/ = 0 and x > 0 are simultaneously satisfied. For y = 0 
and X < 0, although the imaginary part in Eq. (O van¬ 
ishes due to the signum function, its real part becomes 
imaginary, i.e. /3 = i \/fx|, which does not correspond 
to propagation SPP modes. Studying the permittivity 
dependence of x and y in Eq. (|3]) and solving the condi¬ 
tion y = 0 with respect to the dielectric gain part 624 for 
£d ^ —Em, we obtain two exact solutions, i.e., £ 20 , —^ £24 
of the form 



(5) 


The result in Eg . (151) is in agreement with the one de¬ 
rived in Refs. [J, l23| following yet a different derivation 
path. Invoking the physical argument of the SPP wave 
bound to the dielectric-metal interface, we read that only 
£ 2 d is of physical relevance, since £ 2 ^ leads to waves ra¬ 
diating in the transverse towards the interface direction 
[^ . Taking into account the, by definition, positive real 
domain of and the dependence of the latter on the 
metal-dielectric components, we read that the following 
inequality has to be satisfied 

IpmlP > 2eid£2m (6) 

For £d = —£m both x and y diverge exhibiting asymp¬ 
totically the same image, thus 

lim JJ[/3] = lim S[/3] —>• -l-oo (7) 

Sd —^ Sd —^ 

This in turn means that the former complex point does 
not belong to the set domain of the lossless SPP propa¬ 
gation since y ^ 0. 

Solving, on the other hand, the equation x = 0 with 
respect to the dielectric gain £24 for £d ^ —£m, we may 
determine the critical value £c distinguishing the regimes 
of lossless and prohibited SPP propagation, namely 


Ec 



- 1 


( 8 ) 


Equating Eqs. © and ([5|), i.e., y = 0 = x, we obtain 
the condition £i 4 £ 2 m (su — £im) = 0 which reduces to 
Eld = Elm, since eid, e2m > 0. Replacing the former 
value of £id in Eqs. (I5|) and (H]) we obtain T Ec and 
^24 = Ec = E2m for £im < £ 2 m and Elm > E2m, respec¬ 
tively. The former case is obviously a contradiction. The 
latter case corresponds to the singularity point in Eq. 
©, where x,y 7^ 0, thus it is a contradiction as well. In 
other words, x and y do not become zero simultaneously, 
implying that the critical value £c is not an element of 
the domain set of the wave number /3(y = 0) = Pq. This 
is in agreement with the propagation length L in Eq. 
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(ED. Indeed, when y = Q then L tends to infinity, which 
means that the SPP wave number must exhibit a nonzero 
value. An even more interesting point, unveiled from the 
y = Q = X analysis, is the estimation of the /3o behaviour 
when approaching the critical point, £c {y = 0 and 

a; —>■ 0), described by Eq. ([7]) for ,3 ^ Po with eim > E 2 m- 
Physically, the former point corresponds to the wave elec¬ 
trostatic character of zero phase velocity, known in lit¬ 
erature as the surface plasmon mode [^. Including the 
discontinuity at the point Ec, the entire codomain of /Sq 
is described as follows 

{ Real, X > Q 

Imaginary, a; < 0 (9) 

Complex Infinity, a; ^ 0 


w, so that by tuning uj we may control the values of 
e-pj- lying below or above £c- Precisely, for interfaces 
comprised of an active dielectric and a Drude metal [2^ , 
£m is given as 






+ r2 




(u;2-Hr2) 


( 10 ) 


where £h denotes the high frequency permittivity, ojp is 
the plasma frequency and P accounts for metal losses in 
frequency units [^ . 

By virtue of Eq. cni) we can express the epj- and 
Ec in terms of the frequency oj. Moreover, taking Eq. 
© into consideration, we may obtain the SPP resonance 
frequency ujgp [ 1 , 0 , 


In the general case of (3, where y (excluding the point 
Ed = —Em) may take nonzero values as well, we observe 
the following. Eor E 2 d < the signum function in 
Eq. ([3]) is negative since then y < 0 implying SJ[/3] < 
0 => L < 0. This means that the imaginary part of 
/3 accounts for losses and the SPP amplitude decreases 
along the propagation surface. Reversely, for E^d < £ 2 d < 
^ 2 d yielding y > 0 ^ sgn(y) > 0, we have S[/3] > 

0 L > 0. In this case the imaginary part of (3 

accounts for gain and the SPPs amplitude increases along 
the propagation surface. In the special case of y = 0 => 
e 2 d = E^d studied above, the SPP amplitude is constant 
along the propagation surface. This behaviour of the 
signum function fully explains the results observed in Ref. 
regarding the SPP amplitude. 

An interesting feature of the lossless SPP propagation 
case, i.e., for E^m < Ec, in regard to the refractive index n 
is that the latter fulfils the condition n{y = 0 ) = n*{y = 
0 ), since its imaginary part vanishes owing to the signum 
function. This in turn, may be considered as the VT 
symmetry phase condition, where n is spatial indepen¬ 
dent. However the structure is not V'T symmetric in the 
narrow sense, the real value of the supported propagation 
constant along the interface admits time-reversal and ge¬ 
ometrical symmetry. Then, the dielectric gain expression 
E 2 d in Eq. ([5]) can be attributed to the VT symmetry 
property satisfied by the lossless SPP propagation and 
denoted as E^d = Epr- We shall keep this denomination 
in what follows. On the contrary, in the case E^m < £c, 
the VT condition is not satisfied, since the refractive in¬ 
dex is imaginary. Subsequently, the critical gain Ec may 
be regarded as the PT-symmetry breaking point of the 
plasmonic system under scrutiny. 

III. ACTIVE DIELECTRIC - DRUDE METAL 
INTERFACE 

It is quite remarkable that Epj- as well as Ec depend on 
the optical properties of the dielectric and metal, that is, 
Eidy Elm a-nd E 2 m- The metal permittivity in turn may 
generally exhibit a dependence on the angular frequency 


I aA 

V Eld + Eh 

It can be proven that for Drude metals the Epj- is always 
smaller than Eq for frequencies lower than Wgp. Thus, 
according to our theoretical results we anticipate loss¬ 
less SPPs propagation for w < ujgp and prohibited SPPs 
propagation for w > uj sp¬ 
in order to verify our theoretical predictions, we cal¬ 
culate the SPP dispersion relation for an interface con¬ 
sisting of silver with Eh = 1, Wp = 1.367 lO^^Hz, 
r = 1.018 10^‘^Hz and silica glass with Eid = 1.69 and 
E 2 d = Epp. The frequency values are conhned in the 
regime imposed by the inequality in Eq. (I6|). In Eig. 
[T]we plot the real (dotted blue line) and imaginary (or¬ 
ange dashed line) part of the normalized SPP dispersion 
relation /3o/fcp {kp = ujpjc) with respect to the normal¬ 
ized frequency ujjujp\ the yellow dash-dot line shows the 
wave number (/c-number) in the dielectric, and the res¬ 
onance frequency ujgp is represented by the horizontal 
green solid line where the interchange between 5R[/3] and 
5[/3] appears. We observe, indeed, that for w < Wsp the 
imaginary part of (3 vanishes while for w > ujsp the SPPs 
wave number is purely real. Subsequently, in the vicinity 
of w = ujap a phase transition from lossless to prohibited 
SPPs propagation is expected (see Section HVl) . 

Pig. [T] highlights the relation between (3 and the 
metal permittivity, and demonstrates the VT symme¬ 
try breaking point Ec, where the 5R[/3] vanishes. In 
Pig. [7] we consider a variable Eid and record the depen¬ 
dence on it of both the magnitude 5R[/3] in Eq. o and 
Epp in Eq. ([5]) for three different frequencies, namely 
ui = {O.dSwp, O.bwp, O.SSwp}. The former is represented 
by color lines on the complex plane defined by ( ei ^, Epp) 
for X and y axis, respectively. Each color line corresponds 
to a different frequency. Eig. [2] unveils that the more 
dense the dielectric is the higher value of the gain we 
need for having undamped SPPs propagation. In ad¬ 
dition, the 5i[/3] vanishes very suddenly as we increase 
the gain, verifying that at this point the VT symmetry 
breaks and the SPPs propagation becomes prohibited. 
According to the aforementioned figure we can tune the 
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FIG. 1: The SPP dispersion relation /3 under VT symmetry 
{e2d = £vt), i-e., Po, with respect to the frequency uj. 51i[/3] 
and ^[/I] are indicated by dotted blue line and orange dashed 
line respectively. The dash-dot yellow line is referred to the 
wave number of light in the dielectric, whereas the horizontal 
solid green line shows the SPP resonance frequency ujap where 
the interchanging between K[/3] and appears, kp = ojp/c 
is used as normalized unit of wavenumbers and ujp as normal¬ 
ized unit for frequencies as well. 


FIG. 2: The VT gain evT with respect to the real part of the 
dielectric eid is plotted for three different frequencies u). The 
colors on each curve represents the magnitude of the SR[/3] nor¬ 
malized to kp. We observe that the more dense the dielectric 
is the more gain is needed to achieve lossless SPP propaga¬ 
tion. For all frequencies there is a critical gain for which the 
K[/I] transmits suddently from high magnitude (bright col¬ 
ored line) to zero value (black line), showing gain saturation 
and a phase a transition from undapted SPPs to forbidden 
propagation. 


magnitude of the 5R[/3] as well as Sc and e-pj- by choosing 
the appropriate dielectric. 


IV. SIMULATIONS 

In this Section, we verify our theoretical predictions of 
Sections El and uni by solving numerically the full sys¬ 
tem of ME in the frequency domain in a two dimensional 
space (2D) for TM polarization electric and magnetic 
fields. The numerical experiments have been performed 
by virtue of the multi-physics commercial software COM- 
SOL. Precisely, we explore the SPPs propagation length 
L with respect to lo on the interface between two semi¬ 
infinity layers, i.e., an active dielectric and a Drude metal, 
recording the desired phase transition from lossless to 
prohibited SPP propagation. We further demonstrate 
the lossless SPPs propagation, analysing the magnetic 
field intensity along the surface of two known in literature 
conhgurations, the Kretschmann-Raether and the Otto 
conhgurations In our numerical experiments 

the frequency uj is confined in the range [O.Swp, O.TSwp] 
with the integration step Aw = O.Olwp. 

Regarding the active dielectric - Drude metal interface 
described in the previous section, we conduct the near- 
held excitation technique [T^ [H, [2^ to excite SPPs on 
the metallic surface. For this purpose, a circular EM 
source of radius R = 20nm has been located lOOnm 
above the metallic surface acting as a point source, since 
the wavelength A of the EM wave in the silica glass is con¬ 


strained to A >> R [2^. In addition. Perfectly Matched 
Layers (PML) are used as boundary conditions. 

In Figl3]we demonstrate, in a log-linear scale, the prop¬ 
agation length L with respect to w subject VT symmetry 
(blue line and open circle). For the sake of comparison, 
we plot L{uj) for the gainless case (green line and hlled 
circles). The solid lines represent the theoretical predic¬ 
tions obtained by Eq. 0, whereas the circles indicate 
COMSOL results. For the numerical calculations, the 
characteristic propagation length has been estimated by 
the inverse of the slope of the Zog(/)jwhere I is the mag¬ 
netic intensity along the interface [1, [13, H^. The red 
vertical dashed line denotes the SPP resonance frequency 
Wgp, in which the phase transition appears. The graphs 
in Fig. [3]indicate that in the presence of the VT gain, i.e. 
£ 2 m = e-PTi the SPPs may travel for very long, practi¬ 
cally infinite, distances. Approaching the resonance fre¬ 
quency Wsp, L decreases rapidly leading to a steep phase 
transition on the SPPs propagation. The deviations be¬ 
tween theoretical and numerical results in Fig. [3] for fre¬ 
quencies near or greater than Wsp are attributed to the 
fact that in the regime ujgp < lo < ujp, there are quasi¬ 
bound EM modes 2^, where EM waves are evanescent 
along the metal-dielectric interface and radiate perpen¬ 
dicular to this. Consequently, the observed EM field for 
w > oJsp does not correspond to SPPs but belongs to the 
quasi-bound modes. 

So far, the theoretical findings in Sections |TT] and 
El have been successfully confirmed. We further pro- 
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FIG. 3: The SPP characteristic length L is demonstrated 
as function of the frequency uj subject VT symmetry (blue 
line/open circles) and in the case of gainless propagation 
(green line/points). In both cases, the solid lines indicate 
the theoretical prediction, whereas the points show the nu¬ 
merical results obtained by COMSOL simulations. A phase 
transition from lossless (large L) to prohibited propagation 
(small L) occurs in the case of VT symmetry. The red verti¬ 
cal dashed line indicates the resonance frequency ujsp, where 
the phase transition takes place. The theoretical curves show 
deviation from simulations for frequencies near and greater 
that ujsp, because in this regime quasi-bound EM modes ap¬ 
pear. 


ceed investigating the VT symmetry in active plasmonic 
systems. We perform COMSOL simulations based on 
the Total Internal Reflection (TIR) method, applied on 
the Kretschmann-Raether and Otto conhgurations, sepa¬ 
rately. Within the former, a thin metal film is sandwiched 
between two dielectrics with the incident wave hitting the 
denser medium. In Otto configuration, the denser the di¬ 
electric and the metal sandwich a lighter dielectric. In 
both configurations a type of silica glass, with dielectric 
constant Ed = 4, is used as denser passive dielectric, 
whereas for active dielectric (lighter medium) as well as 
for metal, we use the materials described in Section nni 
Furthermore, for the COMSOL simulations [H, [s^, we 
utilize a monochromatic plane wave source of frequency 
/ = 870 THz, amounting to 40% of plasma frequency 
ujp, and corresponding to a wavelength A = 345nm in 
the ultraviolet regime. Again, PML are used as bound¬ 
ary conditions. 

According to the Eq. ([5]), the VT gain for which SPPs 
propagate without losses is calculated to be e-pj- = 0.012. 
We observe, however, that in the numerical experiments 
a larger gain is needed (of the same magnitude though), 
namely ipj- = 0.026, for both configurations. This devia¬ 
tion from the theoretical value can be justified if one takes 
into consideration that the assumption of semi-infinitely 
thick metal and dielectric layers composing the interface 
[ 2 ^, under which the SPP dispersion relation of Eq. (Cl) 


holds true, is experimentally not fully satisfied. 



z 


FIG. 4: COMSOL simulations show the intensity distribution 
of the magnetic field on a Kretschmann-Raether configura¬ 
tion, for (a) £ 2 d = 0 (No gain) (b) e^d = £vt = 0.026 {VT 
symmetry). 


In the Kretschmann-Raether configuration a thin 
metal of thickness d = A5nm has been used for exciting 
SPPs. The resulting propagation is illustrated in Fig. [4^ 
under lack of gain and in Fig. [Hd for the VT symmet¬ 
ric case. The corresponding profiles of the magnetic field 
intensity along the interface are demonstrated in Fig. [51 
where the lossless SPP propagation is evident. 



0 5 10 15 20 25 30 

X (/r m) 


FIG. 5: Characteristic profile of the magnetic field intenstity 
along the interface at Kretschmann-Raether configuration of 
Fig. m for (a) e 2 d = 0 (No gain) (b) S 2 d = £vT = 0.026, 
where lossless SPP propagation is achieved. 

In the Otto configuration, on the other hand, the SPPs 
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excitation is succeeded by means of an active dielec¬ 
tric of thickness d = ISOnm which has been used be¬ 
tween a non-active dielectric and a metal. By analogy to 
the Kretschmann-Raether configuration experiment, we 
present the SPP propagation without gain in Fig. Ha) 
and for the VT gain with e-pj- = 0.026 in Fig. Hb)- In 
Fig. [3 the corresponding profiles of the magnetic inten¬ 
sity along the interface are presented, unveiling again a 
clear lossless SPPs propagation in the VT case. 



z 


FIG. 6: COMSOL results for the intensity distribution of the 
magnetic field at an Otto configuration, for (a) e 2 d = 0 (No 
gain) (b) e^d = £VT = 0.026 {VT symmetry). 




X (/r m) 

FIG. 7: Characteristic profile of the magnetic field intenstity 
along the interface at Otto configuration of Fig. [6l for (a) 
S 2 d = 0 (No gain) (b) e 2 d = ivT = 0.026, where lossless SPP 
propagation is achieved. 


V. CONCLUDING REMARKS 

Summarizing, we have investigated the role of ac¬ 
tive/gain dielectrics in plasmonic systems. In particular, 
we have studied the propagation properties of surface 
plasmon polaritons (SPPs) along an interface confined 
by two semi-infinite layers: a dielectric and a metal. We 
have calculated an exact expression £ 2 ^ for the dielectric 
gain £ 2 d: for which the metal losses have been completely 
counterbalanced, resulting to lossless SPPs propagation 
along the interface. We argued that the a plasmonic sys¬ 
tem characterized by the aforementioned lossless prop¬ 
agation may be related to VT symmetric systems, i.e., 
epj- = £ 2 ^. Within the VT symmetry, a critical gain 
£c exists distinguishing between the real and imaginary 
part of the SPP dispersion relation. This distinction cor¬ 
responds to a phase transition from lossless to prohibited 
SPP propagation. It is remarkable that the £pj- as well 
as the £c depend on the optical properties of the interface. 

We applied our theory to interfaces consisting of Drude 
metals and gain dielectrics demonstrating the predicted 
by the theory lossless propagation as well as the phase 
transition at the SPP resonance frequency Wsp- We per¬ 
formed numerical simulations with COMSOL software, 
using the near-held excitation method in order to inves¬ 
tigate our theory, verifying successfully all the theoretical 
predictions. We also performed COMSOL simulations for 
two different plasmonic conhgurations based on the TIR 
method- Kretschmann-Raether and Otto conhgurations- 
where lossless SPP propagation can be achieved. 

Active metamaterials may be designed to have the 
desirable frequency-dependent permittivity response, as 
Eq. (H points out; these metamaterials could be used 
for the fabrication of VT symmetric plasmonic systems, 
providing inhnite SPPs propagation. The active meta¬ 
materials may be used to design VT symmetric plas¬ 
monic integrated circuits which could transfer informa¬ 
tion in sub-wavelength scales for large (theoretically in¬ 
hnite) distance, rather than a passive plasmonic system 
where SPPs propagate for few micrometers. Moreover, 
we demonstrated that there is a threshold in the VT gain 
values, above which the VT symmetry breaks and thus 
the system passes from lossless to prohibited propaga¬ 
tion. The gain threshold as well as the VT gain depend 
on the optical properties of the dielectric and metal, sub¬ 
sequently we could control the SPPs propagation by tun¬ 
ing the dielectric constant of metal £m or the real part of 
dielectric permittivity ei^; for instance, the former, i.e. 
£m, is usually frequency-dependent, thus we can inter¬ 
change between lossless and prohibited SPP propagation 
by tuning the frequency of the incident EM wave. 
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